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ABSTRACT 

Radio  positions  of  the  three  sources  have  been  determined  with 
the  two  90-ft  antennas  vorking  as  an  interferometer  with  an  rms  ac¬ 
curacy  in  both  coordinates  better  than  10  sec  of  arc.  Direct  photo¬ 
graphs  show  that  a  star-like  object  exists  within  the  error  rectangle 
at  each  of  the  source  positions.  Exceedingly  faint  wisps  of  nebulo¬ 
sity  are  associated  with  the  stars  in  3C  48  and  3C  196.  The  observa¬ 
tions  are  incomplete  for  3C  286  in  this  regard.  Photoelectric  photome¬ 
try  of  the  stars  shows  each  to  have  quite  peculiar  color  indices,  most 
closely  resembling  the  colors  of  old  novae,  or  possibly  white  dwarfs, 
but  we  are  not  suggesting  identification  with  these  types  of  stars. 
Photometry  of  3C  48  through  13  months  shows  the  star  to  be  variable  by 
at  least  AV  «*  0?4.  The  radio  flux  appears  to  be  constant.  Optical 
spectra  for  3C  48  show  several  very  broad  emission  features,  the  most 
intense  at  X3832  being  unidentified.  Spectra  by  Schmidt  of  3C  196  and 
3C  286  show  other  unusual  features.  The  radio  structure  of  the  three 
radio  stars  are  similar  in  that  each  haB  an  unresolved  core  of  <1" 
diameter.  However,  3C  196  and  3C  286  show  halos  of  12"  and  20"  respec- 
j  tively,  while  no  radio  halo  has  been  detected  for  3C  48. 

It  is  shown  that  the  radiant  flux  in  the  optical  region  can  be 
computed  from  the  radio  flux  data  and  the  theory  of  synchrotron  radia¬ 
tion  for  3C  48  and  3C  196,  but  not  for  3C  286.  This,  together  with 
other  arguments,  suggests  that  the  optical  as  well  as  the  radio  flux 
could  be  due  to  the  synchrotron  mechanism,  but  the  arguments  are  not 
conclusive. 

We  have  used  the  assumption  of  minimum  total  energy  to  compute  the 
energy  in  relativistic  particles  and  magnetic  field  required  by  the 
synchrotron  mechanism  to  explain  the  observed  emission.  The  magnetic 
field  in  each  of  the  core  components  is  near  0.1  gauss  and  depends 
mainly  on  the  assumed  angular  size  of  the  emitting  region.  The  total 
energy  in  the  core  components  is  near  10^°  ergs.  The  rate  of  radiation 
is  such  that  the  energy  in  relativistic  electrons  must  be  replaced  in  a 
time  scale  of  a  few  years  if  the  value  of  the  magnetic  field  determined 
in  this  way  is  correct. 

The  frequency  of  occurrence  of  radio  stars  is  examined  and  they 
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are  estimated  to  comprise  from  5  to  10$  of  sources  in  the  3C  catalogue. 
The  percentage  is  likely  to  be  less  for  fainter  sources.  Rough  limits 
have  been  estimated  for  the  mean  distances  of  these  radio  stare.  A 
mean  distance  of  approximately  100  pc  is  suggested. 

I.  INTRODUCTION 

One  of  the  major  programs  of  the  Owens  Valley  Radio  Observatory 
of  the  California  Institute  of  Technology  is  the  determination  of  pre¬ 
cise  positions  of  discrete  radio  sources.  The  radio  observations  are 
made  with  the  two  90-ft  antennas  working  as  an  interferometer  at 
several  spacings  ranging  from  200  ft  to  1600  ft.  The  E-W  direction  is 
used  to  determine  right  ascension,  and  the  N-S  direction  to  find  decli¬ 
nation.  The  observational  technique  for  declination  measurements  has 
been  described  by  Read  (1962)  and  the  entire  problem  and  results  will 
be  discussed  elsewhere  by  Matthews  and  Read.  Errors  in  determination 
of  position  in  both  right  ascension  and  declination  can  now  be  made 
smaller  than  5  seconds  of  arc  under  favorable  conditions.  With  this 
high  positional  accuracy,  the  search  for  optical  identification  is  now 
much  more  efficient  than  similar  searches  made  several  years  ago,  and 
a  number  of  new  identifications  have  already  been  made  (Bolton  I960; 
Maltby,  Matthews,  and  Moffet  1963,  Matthews  and  Schmidt,  unpublished). 

Identifications  to  date  by  all  workers  have  shown  that  radio 
sources  are  associated  with  galactic  nebulae,  supernovae  remnants,  and 
external  galaxies  both  "normal"  and  peculiar.  The  distribution  of  dis¬ 
crete  sources  above  b  =  ±.  20°  is  isotropic  and  has  usually  been  attri¬ 
buted  to  galaxies  alone.  No  star,  except  the  sun,  has  previously  been 
identified  with  a  radio  source.  The  purpose  of  this  paper  is  to  present 
evidence  for  the  identification  of  three  radio  sources  with  objects 
which  are  star-like  in  their  appearance  on  direct  photographs,  and  in 
their  photometric  and  spectroscopic  properties. 

II.  RADIO  AND  OPTICAL  PROPERTIES  OF  THE  THREE  SOURCES 

Our  attention  was  drawn  to  3C  48,  3C  196,  and  3C  286  as  peculiar 
radio  objects  because  of  their  high  radio  surface  brightness.  Measure¬ 
ments  of  the  brightness  distribution  (Maltby  and  Moffet  1962)  along 
both  a  N-S  and  E-W  baseline  at  the  Owens  Valley  Radio  Observatory  with 
a  maximum  baseline  of  1600  wavelengths  showed  that  these  three  sources 
are  single  with  radio  diameters  of  less  than  30  seconds  of  arc.  The 
Jodrell  Bank  observations  of  brightness  distribution  with  four  base¬ 
lines  from  2200X  to  6l,000\  (Allen,  Anderson,  Conway,  Palmer,  Reddish 
and  Row'son  1962)  have  shown  that  even  at  the  longest  baseline  of 
61,0001,  3C  48  Is  unresolved  in  the  E-W  direction,  which  means  that  the 
radio  diameter  is  less  than  1  sec  of  arc  E-W.  Rowson  (1962)  has  also 
shown  that  the  diameter  is  less  than  1  sec  of  arc  in  the  N-S  direction. 
However,  the  Jodrell  Bank  observations  do  show  some  structure  in  3C  196 
and  3C  286  in  the  E-W  direction.  The  simplest  two- component  model 
fitting  the  E-W  intensity  distribution  for  3C  196  is  that  75$  of  the 
flux  comes  from  a  halo  of  about  12"  diameter,  while  the  remaining  25$ 
of  the  flux  is  in  an  unresolved  core  of  less  than  1"  diameter.  For 
3C  286,  40$  of  the  flux  comes  from  a  halo  of  diameter  ^20" ,  and  the  re- 
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malning  6056  ia  again  in  an  unresolved  core  of  diameter  leas  than  1". 

We  are  indebted  to  H.  P.  Palmer  for  the  data  prior  to  publication, 
upon  which  these  diameters  are  based. 

These  small  radio  diameters  together  with  the  large  observed  radio 
flux  initially  suggested  that  the  three  sources  might  be  additional 
examples  of  distant  galaxies  of  large  red  shift  such  as  3C  295,  which 
shows  a  similar  radio  surface  brightness.  Consequently,  when  precise 
radio  positions  were  available,  direct  photographs  were  made  of  each 
field  with  the  200-inch  telescope  in  the  near  red  spectral  region 
(103a  E  plates  plus  Schott  RG1  filter). 

The  first  object  studied  was  30  48  (Matthews,  Bolton,  Greenstein, 
Mtlnch,  and  Sandage  1961).  A  direct  plate  was  taken  on  September  26, 
i960,  with  every  expectation  of  finding  a  distant  cluster  of  galaxies, 
but  measurement  of  the  plate  gave  the  unexpected  result  that  the  only 
object  lying  within  the  error  rectangle  of  the  radio  position  was  one 
which  appeared  to  be  stellar.  The  stellar  object  was  associated  with 
an  exceedingly  faint  wisp  of  nebulosity  running  north-south  (surface 
brightness  ~-23  mag/arcsec^  in  V)  and  measuring  12"  by  5"  (N-S  x  E-W). 

The  stellar  object  lies  about  3"  north  of  the  center  of  the  nebulosity. 
The  peculiarity  of  the  nebulosity,  together  with  the  excellent  agree¬ 
ment  between  the  radio  position  and  the  optical  object,  made  it  almost 
certain  that  an  identification  had  been  achieved.  But  the  nature  of 
the  optical  source  remained  in  doubt  because  in  late  i960  the  existence 
of  radio  stars  was  not  generally  considered  a  serious  possibility. 

Two  spectrograms  were  taken  with  the  prime  focus  spectrograph  at 
the  200-in.  on  October  22,  i960.  One  covered  the  blue-green  region 
from  \3iOO  to  X5000  with  a  dispersion  of  400  A/mm.  The  other  covered 
the  region  from  X3100  to  about  X7OOO  on  an  Eastman  103a  F  plate  with  a 
dispersion  of  800  A/mm.  The  blue-violet  spectrum  is  extremely  peculiar, 
the  only  prominent  features  being  several  strong,  very  broad  emission 
lines.  The  three  strongest  occur  at  X4686  (intensity  4),  X4580  (2)  and 
13832  (6).  The  broad  emission  line  at  X3832  is  the  most  striking 
feature  and  as  yet  has  not  been  identified.  The  most  obvious  identifi¬ 
cation  of  the  X4686  line  is  with  He  II.  If  this  is  correct,  then  the 
measured  wavelength  of  X4686.2  +  1  shows  that  the  radial  velocity  of 
the  object  must  be  less  than  100  km/sec.  The  lines  could  not  be  identi¬ 
fied  with  any  plausible  combination  of  red-shifted  emission  lines.  The 
total  width  of  the  two  strongest  lines  at  half  intensity  points  is 
about  22  A  for  X4686  and  about  30  a  for  X3832.  The  half,  half-widths, 
expressed  in  km/sec,  would  indicate  a  velocity  field  (either  random  or 
systematic)  within  which  the  emission  lines  are  formed  of  about 
1200  km/sec  for  the  X3832  line,  and  J00  km/sec  for  the  X4686  line.  No 
strong  emission  lines  are  present  in  the  red  although  several  faint  ones 
do  exist.  In  particular,  Ha  is  definitely  absent.  Spectrograms  of 
higher  dispersion  were  subsequently  obtained  by  Greenstein  and  Mtinch 
and  a  complete  discussion  of  the  spectroscopic  features  will  be  given 
by  them. 

Photometric  observations  of  the  3C  48  optical  object  confirm  its 
peculiar  nature.  On  October  23,  i960,  the  photometry  gave  V  ■  16.06, 
B-V=  O.38,  U-B  ■  -0.6l,  colors  which  are  similar  but  not  identical  to 


old  novae  (Walker  1957)/  and  to  some  white  dwarfs  (Greenstein  1958), 
but  are  quite  different  from  ordinary  stars  and  galaxies.  This  point 
will  be  discussed  later  in  this  section. 

An  effort  was  made  in  the  case  of  3C  48  to  resolve  the  optical 
image.  On  a  night  of  good  seeing  a  series  of  exposures  ranging  from 
10  minutes  to  15  seconds  was  made  at  the  200-in.  prime  focus  (scale  » 
11.06  arcsec/mm)  on  Eastman  103a  0  plates.  On  the  shortest  exposure 
plate  (15s)  the  image  diameter  of  3C  48  was  measured  to  be  0.09  mm, 
which  corresponds  to  1"  are.  This  is  the  same  diameter  as  images  of 
stars  of  the  same  apparent  brightness  on  the  plate.  The  image  of 
3C  48  on  all  plates  is  sharp  and  appears  to  be  stellar. 

A  second  epoch  Sky  Survey  plate  was  taken  by  W.  C.  Miller  on 
January  18/19,  1961,  with  the  48-in.  Schmidt  to  check  for  a  detectable 
proper  motion.  This  plate  was  centered  identically  with  the  base  plate 
0  30  of  the  original  Sky  Survey  taken  on  December  2l/22,  1949,  giving 
an  11-year  interval.  Inspection  of  the  two  plates  in  a  blink  compara¬ 
tor  showed  no  detectable  proper  motion  relative  to  neighboring  compari¬ 
son  stars.  The  proper  motion  is  less  than  o'.'05/yr  (a  value  which  could 
have  been  detected  by  this  method). 

Optical  photometry  of  3C  48  continued  sporadically  during  1961 
with  the  results  given  in  Table  1.  The  most  striking  feature  of  these 
data  is  that  the  optical  radiation  varies'.  Unfortunately  our  time  re¬ 
solution  is  very  poor.  The  only  evidence  for  short  term  fluctuations 
is  the  data  obtained  on  October  ll/l2,  1961,  where  the  observations 
listed  were  made  in  a  time  interval  of  15  minutes.  We  believe  that 
the  observed  difference  in  intensity  between  these  times  is  probably 
real  because  the  local  standard  star  D  (see  Fig.  2)  was  observed  after 
each  of  the  3C  48  observations.  The  differences  between  the  measure¬ 
ments  of  D  were  only  AV  «  CW007 ,  A(B-V)  «  0?007  and  A(U-B)  -  0?010, 
whereas  the  3C  48  differences  are  AV  =  0?044,  A(b-V)  =  0?053  and 
A(.U-B)  =  0?024. 

The  probable  errors  fbr  most  of  the  data  of  Table  1  are  less  than 
±  0?02  except  for  the  60"  telescope  data  of  November  19/20,  i960,  where 
the  error  is  ±  0?04.  We  therefore  conclude  that  the  night-to-night 
variations  and  those  over  the  13-month  interval  are  undoubtedly  real. 

It  remains  for  future  observations  to  determine  if  short-term  fluctua¬ 
tions  with  a  time  scale  of  minutes  exist  in  the  optical  radiation  simi¬ 
lar  to  those  found  and  studied  in  old  novae  by  M.  Walker  (1954,  1957)* 
The  work  of  Smith  and  Hoffleit  (1961)  shows  that  there  has  been  no 
large  systematic  variation  in  the  magnitude  of  3C  48  over  the  last  60 
years. 


Special  observations  for  the  constancy  of  the  radio  flux  at 
1420  Mc/s  were  made  at  the  Owens  Valley  by  K.  Kellermann  for  more  than 
a  week  in  February  1962,  to  see  if  similar  fluctuations  are  present  in 
the  radio  flux- data.  No  variation  greater  than  3$  was  detected,  the 
probable  error  of  a  single  observation  being  about  1.2$.  Comparison 
with  flux  levels  determined  four  months  earlier  showed  no  variation 
within  the  probable  error  of  about  5$-  Although  the  data  are  not  con¬ 
clusive,  our  present  data  suggest  that  the  radio  flux  of  3C  48  is  con- 
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stant  while  the  optical  flux  varies  by  at  least  a  factor  of  1.4  in  an 
interval  of  13  months. 

Following  the  identification  of  3C  48,  the  other  two  sources  were 
found  in  the  same  way.  For  both  3C  196  and.  3C  286,  the  only  objects  on 
the  direct  photographs  within  the  error  rectangles  have  a  stellar  ap¬ 
pearance.  Furthermore,  the  stellar  object  in  the  position  of  3C  196  is 
also  associated  with  an  exceedingly  faint  wisp  of  nebulosity  which  ap¬ 
pears  as  a  tail  in  position  angle  120°  on  the  stellar  image.  The  wisp 
is  slightly  curved  and  has  dimensions  of  only  3"x  1"  so  it  is  smaller 
than  the  nebulosity  associated  with  3C  48.  The  nebulosity  near  3C  196 
is  extremely  difficult  to  see  on  a  103a  D  plate  behind  a  GG  11  filter 
exposed  for  50  minutes,  but  is  easily  seen  on  a  plate  taken  for  us  by 
W.  A.  Baum  on  103a  E  emulsion  behind  red  plexiglas  for  60  minutes.  We 
have  no  plates  of  3C  48  taken  on  103a  D  emulsion  for  comparison,  but  on 
a  103a  F  plate  behind  Schott  RG2  (126m  exp)  and  on  a  103a  E  +  Schott 
RG2  (90m  exp)  as  well  as  on  several  103a  0  plates,  the  3.C  48  nebulosity 
is  visible  and  stronger  than  that  associated  with  3C  196.  We  have  not 
detected  nebulosity  around  3C  286,  but  again  the  only  available  plate 
is  a  103a  D  +  GG  11  which,  in  view  of  the  extreme  difficulty  of  detec¬ 
ting  the  wisp  in  3C  196  with  this  emulsion  and  filter  combination,  does 
not  allow  us  to  conclude  that  3C  286  has  no  nebulosity. 

The  colors  and  magnitudes  of  the  stellar  objects  at  the  radio  posi¬ 
tions  of  3C  196  and  3C  286  were  measured  photoelectrically  with  the 
200-in.  and  were  found  to  be  quite  unusual,  as  was  3C  48.  The  available 
data,  determined  on  only  one  night  for  each,  are  listed  in  the  last  two 
lines  of  Table  1.  The  colors  of  all  three  sources  are  plotted  in  Fig.l, 
which  also  shows  three  solid  lines:  (C)  the  normal  main  sequence  rela¬ 
tion  for  luminosity  class  V  stars,  (B)  black  body  radiators  and  (A)  ob¬ 
jects  whose  energy  distribution  per  unit  frequency  interval  is  of  the 
form  F(v)oc  v-n.  Lines  are  marked  on  curve  A  for  various  n  values 
ranging  from  n  =  0.0  to  n  =  2.0  in  steps  of  0.2.  This  form  of  the 
energy  spectrum  is  that  predicted  for  synchrotron  radiation  under  cer¬ 
tain  conditions  as  explained  in  Section  III.  The  lines  B  and  C  were 
computed  by  the  methods  of  Appendix  A.  The  black  body  line  differs 
slighti.y  from  that  computed  by  Arp  (1961)  for  reasons  explained  in  Ap¬ 
pendix  A.  Note  that  the  three  radio  sources  fall  close  to  line  A.  The 
only  other  known  stellar  sources  in  this  general  region  of  the  diagram 
are  old  novae  (Walker  1957,  Table  1)  and  a  very  few  white  dwarfs.  Most 
white  dwarfs  fall  on  or  below  the  black  body  line.  Those  old  novae 
which  do  not  have  composite  spectra  are  plotted  as  crosses  in  Fig.  1  to 
illustrate  the  point.  We  are  indebted  to  R.  Kraft  for  unpublished 
spectral  data  on  the  old  novae. 

M.  Schmidt  has  obtained  spectra  of  both  3C  196  and  3C  286  and 
their  spectral  peculiarities  again  confirm  that  we  are  dealing  with  un¬ 
usual  stellar  objects.  Schmidt  reports  that  at  a  dispersion  of 
400  A/mm,  3C  196  has  a  continuous  spectrum  with  no  prominent  features 
either  in  emiss-ion  or  absorption  (contrary  to  3C  48).  He  also  reports 
(1962a)  from  several  spectrograms  that  3C  286  has  a  spectrum  shoving  one 
emission  line  and  possibly  one  absorption  feature.  The  observed  fea¬ 
tures  are  not  the  same  as  those  observed  in  3C  48. 
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Finally,  a  second  epoch  Sky  Survey  plate  of  3C  286  was  obtained 
with  the  48"  Schmidt  on  April  6/7,  1962,  by  J.  Berger.  Comparison  in 
the  blink  comparator  with  the  base  Survey  plate  0  131  taken  on  June  5 /6, 
1950,  showed  that  3C  286  has  no  detectable  proper  motion  over  the  11- 
year  interval.  This  puts  an  upper  limit  to  the  motion  of  0"05/yr. 

Table  2  summarizes  the  radio  and  optical  position  measurements  now 
available  for  these  three  sources.  The  quoted  radio  positions  are  those 
which  were  determined  before  the  optical  identifications  were  made. 

The  radio  right  ascension  measurements  are  due  to  T.  A.  Matthews,  the 
radio  declination  measurements  are  due  to  R.  B.  Read  (1962).  After  the 
preliminary  measurements  of  the  200"  plates  for  optical  position  were 
made  by  Matthews  to  identify  the  stars,  the  final  measurements  of  the 
plates  were  done  by  R.  F.  Griffin  and  the  precise  optical  positions  of 
Table  1  are  due  to  him.  The  agreement  between  the  radio  and  optical 
positions  is  excellent  and  strongly  substantiates  the  reality  of  the 
identifications.  The  quoted  errors  of  Table  2  are  rms  values. 

Figures  2,  3  and  4  are  finding  charts  for  the  optical  objects.  In 
Fig.  1  for  3C  48  the  two  photometric  comparison  stars  B  and  D  are  shown. 
They  have  been  used  as  local  standards  in  the  photometry  done  after 
August  17/18,  1981,  and  have  adopted  values  for  star  B  of  V  ■  13«53» 

B-V  =  0.50,  U-B  =  0.00;  and  for  star  D  of  V  =  14.54,  B-V  »  0.66,  U-B  - 
0.05.  Local  standards  will  be  determined  around  3-C  196  and  3C  286  in 
the  future  when  these  objects  are  followed  more  closely. 

HI.  POWER  SPECTRUM  DATA 
a)  The  Radio  Data 

We  have  compiled  and  tabulated  the  available  radio  flux  data  for 
the  three  sources  in  Tables  3>  4  and  5*  The  measurements  are  from  many 
sources,  either  in  the  literature,  or  from  unpublished  data  from  the 
Owens  Valley  Radio  Observatory  generously  made  available  by  Kellermann 
and  Bartlett.  The  references  are  listed  at  the  bottom  of  Table  5>  The 
flux  values  are  in  units  of  10 ■^“watts (meter ) ( c/s ) and  have  been 
taken  directly  from  the  various  authors  with  no  attempt  to  put  the 
various  absolute  calibrations  on  a  strictly  homogeneous  system.  The 
power  spectra  in  the  radio  range  from  160  Me  to  3200  Me  are  shown  in 
Figures  5i  6  and  7*  Probable  errors  are  given  in  the  tables  where 
known.  As  a  guide,  the  size  of  a  10$  probable  error  is  shown  in 
Figures  5,  6  and  7. 

The  data  are  most  complete  for  3C  48,  and  it  is  evident  from  Fig. 

5  that  no  single  straight  line  can  be  drawn  through  all  the  points. 

The  best  compromise  single  spectral  index  which  satisfies  all  of  the 
points  except  the  three  in  the  3000  Me  region  is  -0.59-  Since  these 
three  high  frequency  observations  carry  so  much  combined  weight,  we 
conclude  that  a  single  spectral  index  is  not  acceptable.  We  have  arbi¬ 
trarily  fitted- segments  of  two  power  laws  to  the„data  with  indexes 
nq  =  -0.47  for  frequencies  lower  than  v  =  7  x  10°  c/s  and  <*  -O.67  for 
v  >7  x  10®  c/s.  Kellermann,  Long,  Allen,  and  Moran  (1962),  using  un¬ 
published  low  frequency  flux  measurements,  believe  that  the  low  fre¬ 
quency  spectral  index  is  even  flatter  than  n  ■  -0.47.  They  have  sug- 
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gested  that  a  high  energy  cut-off  is  responsible  for  the  curvature  seen 
in  the  radio  spectrum.  However,  in  view  of  the  optical  data  to  be  dis¬ 
cussed  later,  we  wish  to  suggest  that  it  might  be  interpreted  as  the 
low  energy  cut-off  of  the  electrons  producing  the  synchrotron  radiation. 
The  possibility  of  a  low  frequency  cut-off  is  not  discussed  further  here 
because  of  the  limited  observational  data. 

The  index  -O.67  fits  the  3C  48  radio  data  for  v  >•  7  x  10®  c/s  quite 
well  and  we  shall  see  later  in  this  section  that  an  extension  of  F(v) 
into  the  optical  region  using  the  synchrotron  radiation  theory  with 
n  =  -O.67  predicts  the  observed  optical  flux  with  remarkable  accuracy. 

Fig.  6  shows  that  the  data  for  3C  196  also  cannot  be  well  repre¬ 
sented  by  a  single  spectral  index,  although  most  of  the  points  (except 
the  three  near  v  =  3000  Mc/s)  do  fall  near  the  dashed  line  of  index 
-O.67.  We  shall  show  later  that  the  optical  data,  if  all*of  the  ob¬ 
served  optical  flux  is  assumed  to  be  due  to  synchrotron  radiation,  re¬ 
quire  n  =  -0.74,  and  this  is  shown  as  a  solid  line  in  Fig.  6  for  v>  4 
x  10°  c/s.  The  behavior  of  the  flux  at  low  frequencies  is  not  well  de¬ 
termined;  one  suggestion  is  shown  by  the  dash-dot  line. 

The  few  data  available  for  3C  286  suggest  that  the  radio  spectral 
index  has  the  very  low  value  of  n  *  -0.28,  but  more  observations  at  low 
and  high  frequencies  are  needed  to  confirm  the  result.  Recent  flux 
data  below  400  Mc/s  measured  at  Cambridge  suggest  that  the  spectrum  of 
3C  286  also  flattens  at  the  low  frequency  end  (Conway,  Kellermann  and 
Long  1962). 

b )  The  Optical  Data 

The  U,  B,  V  values  of  Table  1,  Section  II,  have  been  converted  to 
absolute  flux  units  by  the  methods  of  Appendix  B  and  have  then  been 
compared  with  the  radio  data. •  For  3C  48  we  have  adopted  mean  values 
V  «  16.20,  B-V  =  0.40,  U-B  =  -0.59  for  this  comparison,  rather  than 
treating  each  of  the  observations  separately.  We  were  fortunate  in 
asking  W,  A.  Baum  to  measure  UBV  magnitudes  of  3C  48  during  his  Palomar 
run  in  December  1961,  because  in  addition  to  these  three  colors,  he 
measured  the  star  in  his  eight-color  syBtem  which  extends  from  3700  fi 
to  9800  8  and  he  has  made  these  data  available  to  us  for  discussion. 

Table  6  summarizes  the  available  optical  data.  The  effective  fre¬ 
quencies  for  the  UBV  observations  are  directly  computed  from  the  effec¬ 
tive  wavelengths  worked  out  in  Appendix  B  and  given  in  Columns  4,  5  and 
6  of  Table  B1  for  different  values  of  the  effective  index  n.  The  value 
of  n  was  estimated  for  each  of  the  three  sources  from  their  observed 
B-V  values  by  comparing  them  with  Column  3  of  Table  Bl.  Baum's  data 
were  converted  to  absolute  units  by  force  fitting  them,  over  the  fre¬ 
quency  interval  in  common,  to  the  absolute  calibration  of  the  U,  B,  V 
points  given  in  the  first  three  rows  of  Table  6.  The  flux  unit  in 
Table  6  is  10“s9 w  m"2( c/s ) "1 . 

The  radio  and  optical  data  are  plotted  together  in  Fig.  8.  The 
straight  line  with  its  family  of  termination  curves  A,  B,  and  C  is, 
under  certain  conditions  (see  below),  the  theoretical  power  spectrum 


for  an  assemblage  of  relativistic  electrons  whose  number  In  energy 
range  E  to  E  +  dE  Is 


N(E)  dE  =  K  E-^  dE 


(1) 


for  E^«  E  cEg  where  Ej_  and  Eg  are  the  low  and  high  energy  cut-offs. 


The  well  known  theory  of  synchrotron  radiation,  summarized  by  Oort  and 
Walraven  (1956),  by  Woltjer  (1958),  by  Burbidge  (1956)  and  others,  pre¬ 
dicts  the  power  at  frequency  v  radiated  by  the 
be  of  the  form  • 
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in  units  of  ergs  sec-1  cm-2  (c/s)"1.  Hj_  is  the  component  of  the  magne¬ 
tic  field  strength  perpendicular  to  the  velocity  vector  of  the  electron; 
L  =  I.608  x  10l3  if  the  electron  energy  E  is  expressed  in  Bev,  L  ■  6.269 
x  10l®  if  E  is  in  ergs,  and  a  relates  the  frequency  of  observation  v 
with  the  low  and  high  critical  frequencies  vq  and  vg  by 
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where  the  critical  frequencies  are  given  by 
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The  function  F(a)  has  been  tabulated  by  Vladimirsky  (19^8),  by  Oort  and 
Walraven  (1956),  and  by  Westfold  (1958).  In  the  following  calculations 
we  shall  assume  that  the  low  frequency  cut-off  vq  is  below  10?  c/s.  To 
compute  the  theoretical  flux  for  frequencies  greater  than  10°  c/s,  we 
can  put  ap  =  00  with  good  accuracy.  In  what  follows  we  also  put  Hj_ 
equal  to  a  constant  and  thereby  ignore  the  possibility  of  magnetic  field 
gradients.  Such  gradients  will  give  a  different  F(v)  distribution  from 
the  constant  H  model,  as  will  a  model  such  as  has  been  postulated  for 
the  Crab  (Oort  and  Walraven  1956),  where  different  regions  of  the  source 
contain  electrons  of  different  energies.  Under  these  assumptions  equa¬ 
tion  (2)  becomes 

.  F(v)  -  3) 
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where  D  is  a  constant  and  where 
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x f  dt •  (6) 

If  the  frequency  of  observation  v  is  far  enough  removed  from  the  high 
frequency  cut-off  (i.e.  ag*0.01  or  V2 ^ 100  v),  then  the  theoretical 
energy  spectrum  will  be  a  power  law  of  index  {/S  - 1)/2,  in  good  general 
agreement  with  Figures  5>  6  and  ?• 

If  we  now  require  that  all  the  optical  flux  of  30  4-8  is  synchro¬ 
tron  emission  alone,  then  Fig.  8  shows  that  the  high  frequency  cut-off, 
vg,  is  near  10^5  c/s  since  the  optical  points  are  fainter  than  the  ex¬ 
trapolated  v-O.67  relation.  l(ag)  has  been  found  by  numerical  integra¬ 
tion  of  equation  (6)  with  £  =  2.43  corresponding  to  a  spectral  index 
-0.67. 

Once  the  value  of  /i  is  fixed,  the  only  free  parameter  in  equation 
(5)  is  the  lower  limit  of  the  integral  a%.  Curves  A,  B  and  C  of  Fig.  8 
are  predicted  spectra  for  critical  frequencies  vg  of  3  x  lol1*  c/s, 

6  x  lO11*-  c/s  and  1015  c/s  respectively.  Fig.  9  shows  more  detailed 
calculations  in  the  optical  range  compared  with  the  data  of  Table  6. 

The  six  curves  are  for  cut-off  frequencies  ranging  from  4  x  10^^  c/s  to 
10i5  c/s.  A  frequency  of  Vg  =  7  *  10^  c/s  provides  a  remarkably  good 
fit  to  all  of  the  data. 

The  detailed  calculations  show  that  deviations  from  the  simple 
power  law  F(v)<xv-n  occur  if  the  frequency  of  observation  is  closer  to 
the  critical  frequency  vg  than  02  =  0.01.  For  02  >0.01  the  [log,  log] 
plot  is  no  longer  a  straight  line,  but  over  small  intervals  in  fre¬ 
quency  F(v)  can  still  be  adequately  represented  by  an  "effective  index 
n"  which  will  be  larger  than  (fi- l)/2.  Over  the  optical  range  of  the  U, 
B,  V  filters  (8.4  x  10i^»vs>5«5  x  lOi^1)  the  effective  index  for  3C  48 
is  about  1.5  (rather  than  0.67  which  holds  for  ag-cO.Ol)  and  this  ex¬ 
plains  why  the  plotted  point  for  3C  48  in  Fig.  1  lies  close  to  and  be¬ 
low  curve  A  near  the  value  for  n  =  1.3.  It  is  important  to  note  that 

although  curve  A  of  Fig.  1  was  computed  using  the  theoretical  synchro¬ 
tron  power  law  of  equation  5  with  ag  =  0,  objects  which  lie  close  to 
curve  A  are  not  necessarily  synchrotron  radiators.  The  close  fit  of 
observations  to  curve  A  only  means  that  over  the  very  small  frequency 
range  between  the  U  and  the  V  filters  the  energy  distribution  can  be 
approximated  by  F(v)ocv-n  with  the  n  appropriate  to  the  particular 
position  on  curve  A. 

Finally,  with  regard  to  Figures  8  and  9  it  should  be  mentioned 
that  the  value  of  the  spectral  index  (j8~ l)/2  is  quite  critical  for  a 
good  fit  of  the  optical  data.  The  first  calculations  for  3C  48  were 
made  with  an  index  of  0.59  as  suggested  by  the  best  overall  fit  to  the 

radio  data  of  Fig.  5.  The  fit  of  those  calculations  to  the  optical 

data  was  much  poorer  than  in  Fig.  9  with  deviations  over  the  U  to  V 
interval  ranging  to  0?3«  Consequently,  in  making  the  fit  of  Fig.  9 
there  are  in  effect  the  two  adjustable  parameters  js  and  vg.  However, 
jS  can  only  be  adjusted  within  the  limits  Imposed  by  the  detailed  radio 
data  of  Fig.  5  and  Table  3* 
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Similar  calculations  vers  made  for  3C  196  and  3C  286  and  are 
shown  in  Figures  10  and  11.  In  the  case  of  3C  196,  Fig.  10  shows  that 
equation  (5)  gives  a  satisfactory  tie-in  between  the  radio  and  optical 
data  if  the  high  frequency  cut-off  is  about  4  x  lO1^  c/s.  But  again  it 
should  be  emphasized  that  this  does  not  in  itself  prove  that  the  opti¬ 
cal  radiation  is  due  to  synchrotron  emission.  However,  Figures  8  and 
10  are  so  similar  that  it  would  be  surprising  if  the  optical  emission 
were  not  somehow  related  to  the  radio  flux.  On  the  other  hand,  Fig.  11 
shows  that  the  situation  must  be  more  complicated  for  3C  286.  We  aan 
derive  a  curve  such  as  B  which  will  pass  through  one  of  the  optical 
points  (which,  of  course,  is  always  possible  as  long  as  the  optical 
flux  is  less  than  or  equal  to  the  v-(^-l)/2  relation  extrapolated  from 
the  radio  region)  but  no  change  of  parameters  can  predict  the  observed 
color  of  3C  286  as  we  have  done  for  3C  48  and  3C  196.  Consequently, 
we  conclude  that  the  optical  radiation  of  3C  286  is  probably  not  due  to 
synchrotron  emission.  The  high  frequency  cut-off  must  occur  at  fre¬ 
quencies  between  10^  and  10^3  e/s  and  most  if  not  all  of  the  observed 
radiation  from  the  3C  286  star  may  be  thermal. 


In  addition  to  the  agreement  of  theory  and  observation  in  Figures 
8,  9,  10  and  11,  an  indirect  argument  might  be  used  to  suggest  that  the 
optical  flux  is  synchrotron  radiation.  The  argument  concerns  the  ex¬ 
planation  of  the  fluctuations  of  the  optical  intensity  of  3C  48.  The 
time  for  an  electron  to  lose  one  half  of  its  original  energy  E  (in  Bev) 
is  given  (Oort  and  Walraven  1956)  by 


tl/2 


days. 

h|  e 


(T) 


The  energy  of  the  electrons  which  emit  the  large  majority  of  the  opti 
cal  radiation  is  related  to  the  magnetic  field  through  equation  (4), 
because  individual  electrons  radiate  mainly  near  their  own  critical 
frequency.  Putting  v  =  7  x  lQl^  in  equation  (4)  gives  E  ■  6.6  H-1/ 

’  Btituting  in  equation  ("' 
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372 


(7)  we  have 


Bev 


days  (for  v  «  7  x  lO1^  c/s) 


(8) 


H 


for  the  decay  time  of  electrons  responsible  for  optical  radiation.  In 
the  radio  frequencies,  the  decay  time  is  very  much  longer.  If  v  * 

1000  Mc/a,  then  equation  (4)  requires  E  =  7*9  x  10*2  Hj1/2  Bev  which, 
from  equation  (7)  gives 
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387 
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days  (for  v 


lo9  c/s). 


(9) 


If  the  synchrotron  radiation  is  steady,  new  electrons  must  be  in¬ 
jected  at  or  accelerated  to  energy  E  at  Just  the  rate  necessary  to  re¬ 
place  electrons  degraded  in  energy  from  E  due  to  their  own  radiation. 

If  this  steady  state  is  not  exactly  maintained,  then  the  distribution 
function  of  equation  (l)  is  not  exactly  maintained,  and  the  radiated 
power  will  fluctuate  with  the  time  scales  of  equation  (7)«  Equations  (8) 
&(9)  show  that  the  corresponding  time  scale  of  variations  for  the  radio 
frequencies  is  840  times  longer  than  the  optical  time  scale,  which  might 
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then  explain  why  the  radio  flux  is  observed  to  remain  constant  over  the 
same  time  intervals  in  which  the  optical  radiation  fluctuates. 

The  amplitude  of  the  fluctuations  will  be  determined  by  the  fluc- 
tuation  in  the  number  of  electrons  radiating,  or  by  any  change  in  the 
mean  magnetic  field  seen  by  these  electrons.  The  number  of  electrons 
being  observed  by  a  radio  receiver  is  900  times  the  number  observed  by 
•  the  V  filter.  This  takes  into  account  the  electron  spectrum  (equation 

l),  the  bandwidths  used  at  the  two  frequencies  (radio,  4v=  2x  10^o/s; 
optical,  A v = 8 . 6 x 10l3c/s )  and  the  bandwidth  into*which  an  electron 
radiates  at  the  two  frequencies.  About  10$  of  the  radio  electrons  comes 
from  electrons  which  previously  radiated  in  the  optical  region,  thus  the 
presence  of  radio  intensity  fluctuations  depends  primarily  on  the  con¬ 
stancy  of  the  injection  mechanisms  at  energies  from  0.1  to  1  Bev. 
Equation  (8)  shows  that  the  optical  flux  variation  on  a  time  scale  of  a 
day  would  require  to  be  about  1  gauss  and  therefore  the  optical  elec¬ 
trons  would  have  an  energy  about  6  Bev  according  to  equation  (4). 

With  the  synchrotron  mechanism  in  mind  we  requested  W.  A.  Hiltner 
to  measure  the  optical  polarization  of  3C  48,  which  he  did  at  the 
McDonald  Observatory  in  late  i960,  and  he  reports  that  any  polarization 
is  less  than  2$.  This  negative  result  may  mean  either  that  the  optical 
radiation  is  thermal  or  that  the  magnetic  fields  are  sufficiently 
tangled  so  that  the  net  polarization  observed  when  averaged  over  the 
region  of  emission  is  immeasurably  small. 

IV.  SOME  PHYSICAL  PARAMETERS  EVALUATED 
PROM  THE  SYNCHROTRON  THEORY 

Further  insight  into  the  physical  conditions  that  may  exist  in 
these  stars  can  be  gained  by  using  the  synchrotron  mechanism  to  calcu¬ 
late  the  magnetic  field  and  minimum  total  energy  required  to  produce 
the  observed  radiation.  In  doing  this,  we  follow  the  computational 
methods  of  Burbidge  (1956}  see  also  Burbidge  and  Burbidge  1957)* 

Parallel  calculations  have  been  carried  out  assuming  first  that  the 
synchrotron  spectrum  extends  into  the  optical  region  as  in  Figures  8,  9, 
and  10,  and  second  that  the  optical  radiation  is  thermal  and  the  high 
frequency  cut-off  occurs  at  v2  =  1010  c/s.  The  total  flux  emitted  by 
the  source  is  given  by 

L  =  4n  r2  f  F(v)  dv  (10) 

'i 

where  r  is  the  distance  (L  is  not  to  be  confused  with  the  constant  in 
equations  2  and  4).  In  making  the  integration  the  flux  spectral  distri¬ 
bution  has  been  assumed  to  be  of  the  form  F(v)ocv_n  between  the  lower 
cut-off  frequency  and  the  upper  cut-off  frequency  V2.  Using  the  dis¬ 
tribution  of  electron  energies  given  by  equation  (l)  and  noting  that  the 
rate  of  radiation  of  a  single  electron  is  dE/dt  =  -2.368  x  10"3  H2  e2 
ergs  sec-1,  then 

/-E  g  _ 

L  -  2.368  x  10-3  H2  K  J  E2  N(E)  dE 

E1 


(11) 
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where  E  is  in  ergs,  L  in  ergs  (sec)*1,  in  gauss. 

The  total  energy  of  relativistic  electrons  is  Ee  and  is  given  by 


We  also  make  the  assumption  (see  Burbidge  1959)  that  the  heavy  particle 
energy,  Ep,  is  100  times  the  electron  energy  Ee. 

The  energy  in  the  magnetic  field 

E„  -  |i  V  (13) 

where  V  is  the  volume  of  the  emitting  region.  If  we  now  assume  that 
the  total  energy  (Em  +  Ee  +  Ep)  is  a  minimum,  then  it  can  be  shown  that 

Em  =  |  <ep  +  Ee )  •  (1*0 

From  equations  (1),  (ll),  (12),  (13)  and  (l4)  we  can  solve  for  the  mini 
mum  total  energy  ET  =  Em  +  Ep  +  Ee  with  the  result  that 

Et  (min.)  =  4.586  x  lO1?  [c  L(l)]4/?  s9/T  r17/T  (15) 

where  S  is  the  angular  diameter  in  arcseconds,  and  where  L(l)  is  the 
emitted  flux  if  the  source  were  at  a  distance  of  one  parsec.  The  value 
of  C  is  a  slow  function  of  spectral  index  n,  and  of  the  upper  and  lower 
energy  cut-offs,  and  is  given  by 


C- 


"  $)' 
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The  magnetic  field  H  which  produces  the  minimum  value  of  ET  is 

H(Et  min.)  =  3-8l9x  10'1*  [c  LU)]2/7  r-2/^  gauss.'  (l6) 


An  interesting  result  of  equation  (16)  is  that  the  magnetic  field  de¬ 
pends  almost  entirely  on  the  angular  size  and  only  slightly  on  the  dis¬ 
tance  . 

The  ratio  T =.  is  a  measure  of  the  time  scale  of  replenishment 

of  the  relativistic  electrons,  or  of  the  lifetime  of  the  source  itself 
if  no  replenishment  occurs.  Under  the  assumptions  that  have  been  made, 
the  total  energy  E>j  is  236  times  Ee. 

Since  we  do  not  have  any  good  estimate  of  the  distance  to  these 


sources  the  values  of  L(l),  V(l),  H(l),  Eij>(l)  and  2”(l)  are  given  in 
Table  7  for  a  distance  of  one  parsec.  The  value  of  the  absolute  visual 
magnitude  Mv  is  also  given  where  the  lower  limit  corresponds  to  assu¬ 
ming  that  all  the  optical  radiation  is  thermal.  Table  8  gives  the 
various  scaling  factors  to  calculate  the  parameters  for  other  distances. 
The  angular  size  of  3C  48  has  only  an  upper  limit  as  do  the  cores  of 
3C  19*5  and  3C  286.  In  order  to  indicate  the  effects  of  size,  the  tables 
have  been  computed  for  two  assumed  sizes  of  1.0  and  0.1  arcsecond. 

The  most  interesting  feature  of  Table  7  is  the  large  value  of  mag¬ 
netic  field  necessary  in  each  of  the  radio  core  components.  The  value 
is  almost  independent  of  the  assumed  cut-off  frequency,  but  does  depend 
on  the  assumed  size.  The  value  of  H  =  0.20  S""'1  r-®/'  gauss  for  3C  48 
is  not  far  from  the  estimate  of  1  gauss  found  necessary  in  Section  IIIc 
to  explain  the  observed  optical  intensity  fluctuations;  this  perhaps 
adds  further  weight  to  the  suggestions  of  optical  synchrotron  radiation. 


For  a  mean  distance  of  100  pc,  as  is  suggested  .in  Section  VI,  for 
3C  48  the  values  of  total  energy  Ej  =  1.1  x  10 s9/f  erg,  and  of  the 
magnetic  field  H  =  0.054  S*°/7  gauss,  are  almost  completely  independent 
of  cut-off  frequency.  The  cores  of  the  three  objects  all  have  similar 
values  for  these  quantities.  The  time  scale  T  depends  on  cut-off  fre¬ 
quency  and  is  r(v2  =  10-*-®'  c/s)  =  86  s9/7  years  or  r(v2  =  7  x  10^  c/s) 

=  2.4  s9/7  years.  Since  S  is  less  than  1",  both  of  these  time  scales 
are  very  short  and  indicate  either  that  the  star  must  be  continuously 
supplying  relativistic  electrons,  or  that  the  production  of  secondary 
electrons  and  the  acceleration  of  particles  are  efficient  processes. 

In  either  case  a  continuing  supply  of  relativistic  electrons  must  be 
made  available  on  a  short  time  scale  by  some  process  if  the  magnetic 
field  we  have  assumed  is  correct.  But  it  should  be  pointed  out  that  the 
requirement  for  replenishment  need  not  hold  if  we  discard  the  require¬ 
ment  of  minimum  total  energy.  If  the  magnetic  field  is  arbitrarily 
made  very  much  smaller  than  computed  above,  t~  can  be  made  arbitrarily 
long  because  r"ocH-3/2.  Thus,  if  later  investigations  prove  it  diffi¬ 
cult  to  find  a  continuing  supply  of  electrons,  T can  be  made  as  long 
as  the  lifetime  of  the  source  by  lowering  H.  Further  investigation  of 
the  optical  variations  might  show  whether  or  not  the  optical  radiation 
is  indeed  due  to  the  synchrotron  mechanism,  and  thus  determine  the 
value  of  H. 


Table  9  gives  the  values  of  the  physical  parameters  discussed  here 
that  might  be  expected  for  an  average  object  at  the  suggested  mean  dis¬ 
tance  of  100  parsecs.  The  mean  absolute  visual  magnitude  would  be 
Mv  =  +12. 

V.  AN  ESTIMATE  OF  THE  PERCENTAGE  OF  RADIO 
STARS  IN  THE  VARIOUS  CATALOGUES 

a )  The  Strongest  Radio  Sources 

In  the  3 C~,  Mills,  and  the  Caltech  catalogues  there  are  24  radio 
sources  north  of  5  =  -47°  whose  flux  at  958  Mc/s  is  greater  than 
19  x  10“”°  W  m“2(c/s)”1,  a  limit  which  is  just  below  the  faintest  of  the 
radio  stars.  This  listing  omits  all  sources  which  are  identified  with 


galactic  nebulosities  (i.e.  H  II  regions,  aupernovse  remnants,  etc. ), 
or  which  probably  belong  to  that  class  (i.e.  low  galactic  latitude, 
radio  diameter  5 3 ' ,  or  thermal  spectrum).  Of  these  24  sources,  13 
have  been  identified  with  galaxies  and  3  are  now  identified  with  stel¬ 
lar  objects,  leaving  8  unidentified  sources. 

Two  of  the  8  are  probably  extragalactic  because  they  have  the 
characteristic  two-component  structure  (Malt-by  and  Moffet  1962)  found 
for  three  quarters  of  all  extragalactic  radio  sources.  An  additional 
two  are  also  probably  extragalactic  for  similar  reasons,  but  the  data 
are  poorer.  This  leaves  4  possible  candidates,  which  together  with 
the  3  radio  stars  already  known,  suggests  that  as  many  as  7  out  of  the 
24  strongest  radio  sources  (excluding  galactic  nebulosities)  might  be 
radio  stars.  Therefore,  at  least  3/24  (or  12$),  but  less  than  7/24 
(or  29$)  of  the  strongest  sources  which  are  not  galactic  nebulosities 
could  be  radio  stars. 


b )  In  the  3C  Catalogue 

The  diameter  measurements  made  by  the  Jodrell  Bank  workers  (Allen 
et  al«  1962)  of  324  sources  in  the  3C  catalogue  provides  another  possi¬ 
ble  means  of  separating  out  more  candidates.  The  diameter  characteris¬ 
tics  of  the  three  known  radio  stars  are  unique  in  that  they  have  a  halo 
(probably  not  much  larger  than  30’)  and  a  considerable  fraction  of  the 
radiation  concentrated  in  an  unresolved  core.  Including  3C  48,  3c  196, 
and  3C  286,  there  are  13  radio  sources  out  of  the  sample  of  324  which 
have  the  necessary  diameter  characteristics.  None  of  these  have  been 
identified  with  galaxies.  (Four  other  sources  might  also  be  included, 
but  the  information  is  not  complete  enough  to  decide  if  they  fulfill 
all  the  criteria.)  These  13  (or  possibly  17 )  sources  suggest  that 
about  4$  (or  possibly  5$)  of  that  portion  of  the  3C  catalogue  which  has 
been  surveyed  for  angular  diameter  could  be  radio  stars.  This  percen¬ 
tage  is  a  lower  limit  because  it  depends  on  the  degree  of  completeness 
of  the  3C  catalogue  at  low  galactic  latitudes  where  the  effects  of  con¬ 
fusion  and  uneven  background  variations  are  strongest.  Of  course  we 
are  certain  that  the  absolute  lower  limit  to  the  percentage  of  radio 
stars  in  the  catalogue  is  1$,  obtained  from  the  three  known  stars  out 
of  a  sample  of  324. 

The  question  now  arises  whether  these  estimates  of  the  percentage 
of  radio  stars  made  above  are  in  contradiction  with  the  findings  of 
the  Cambridge  group  (Scott,  Ryle,  and  Hewish  1961}  Ryle  and  Clarke 
1961)  that  less  than  1$  of  the  sources  brighter  than  -'10-27  W  m-2 
(c/s)--1-  at  178  Mc/s  can  be  halo  objects  in  the  Galaxy.  The  special 
Cambridge  observations  were  made  in  two  selected  regions  toward  and 
away  from  the  hemisphere  containing  the  galactic  center  (Scott  et  al. 
used  -'120°,  bII~-54°  and  -'300°,  bll  — +63 Ryle  and  Clarke 

used  146°,  b-^-'-jo0  and  JZ  **  -^330°»  b^/v+590  and  separated  by 

about  l8o°  on  the  plane  of  the  sky.  If  an  appreciable  number  of  the 
faint  sources  -detected  in  that  study  had  been  in  the  galactic  halo, 
then  a  marked  difference  in  the  number  of  sources  in  these  two  direc¬ 
tions  would  have  been  found  because  of  the  different  volumes  of  the 
galactic  halo  surveyed.  No  difference  was  observed  greater  than  the 
statistical  fluctuations  of  10+  sources  per  steradian  to  the  limit  of 
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10*27  w  m"2  ( c/s ) "I .  Therefore,  less  than  1$  of  the  radio  sources  to 
this  limit  can  belong  to  the  galactic  halo.  These  findings  are  not  in 
contradiction  with  our  estimates  for  the  following  reasons.  (l)  Our 
estimate  of  4^  to  5 %  as  radio  stars  refers  to  sources  brighter  than  the 
limiting  flux  of  the  3C  catalogue,  8  x  10”26  w  m~2  (c/s)*l.  It  is  very- 
likely  that  this  percentage  will  drop  in  catalogues  which  survey  to 
fainter  flux  levels  because  the  background  count  of  external  sources 
will  begin  to  outnumber  the  galactic  sources  in  the  same  way  that 
galaxies  outnumber  stars  as  we  photograph  to  faint  levels  in  optical 
astronomy.  (2)  The  Cambridge  conclusion  refers  to  objects  in  the 
galactic  halo,  whereas  the  radio  stars  may  be  concentrated  toward  the 
galactic  plane  and  no  evidence  of  anisotropy  would  then  be  detected 
by  their  experiment. 


VI.  ESTIMATE  OF  THE  DISTANCE  OF 
THE  RADIO  STARS 

Two  arguments  may  be  used  to  3how  that  the  three  radio  stars  iden¬ 
tified  so  far  are  probably  nearby. 

a)  The  thirteen  possible  radio  stars  discussed  in  Section  V  have 
a  random  distribution  in  both  latitude  and  longitude  with  no  concentra¬ 
tion  toward  one  hemisphere  of  the  sky.  Assuming  that  these  13  candi¬ 
dates  are  actually  radio  stars,  then  their  roughly  isotropic  distribu¬ 
tion  shows  that  they  cannot  be  distant  halo  objects  because,  if  they 
were,  a  marked  asymmetry  toward  the  hemisphere  containing  the  galactic 
center  would  occur  due  to  the  sun's  peripheral  position  (note  e.g.  the 
globular  clusters).  Therefore  the  radio  stars  mu3t  be  nearby  and  in 
the  galactic  plane. 

b)  A  rough  upper  limit  to  the  average  distance  between  radio  stars 
can  be  estimated  as  follows.  We  3hall  assume  that  the  number  of  radio 
stars  in  the  3C  catalogue  is  0.07  x  471  =  33  sources.  This  number  must 
be. corrected  for  the  portion  of  the  sky  unobserved  in  the  catalogue 
which,  if  the  radio  stars  are  concentrated  toward  the  plane,  gives  an 
incompleteness  factor  of  1.6.  There  would  then  be  about  50  radio  stars 
in  a  whole  sky  survey  carried  to  a  flux  density  greater  than  6  x  10-26 

W  m-2  (c/s)-l  at  178  Mc/s.  An  extreme  assumption  which  leads  to  a 
generous  upper  limit  to  the  average  separation  is  that  the  entire  volume 
of  the  galaxy  is  surveyed  for  these  stars  to  this  flux  level.  Assuming 
a  volume  of  10H  pc3  for  the  galaxy  gives  an  average  density  of  one 
radio  star  per  2  x  109  pc3.  If  we  divide  the  galaxy  into  cubical  cell3 
of  length  d  pc  on  a  side  such  that  each  cell  contains  one  radio  star  at 
its  center','  then  d3  =  2  x  109  pe3  or  d  =  1300  pc.  This  is  of  the  order  a 
of  the  average  distance  of  any  given  radio  star  from  its  nearest  neigh¬ 
bor.  If  the  sun  is  equidistant  between  its  two  nearest  radio  stars, 
then  the  average  distance  between  us  and  the  nearest  radio  star  is 
d/2  =  650  pc. 

Thi3  is  an  extreme  upper  limit  because  we  have  seen  that  the  13 
possible  candidates  for  radio  stars  do  not  have  an  apparent  concentra¬ 
tion  toward  the  galactic  plane  and  have  no  asymmetry  in  their  distribu¬ 
tion  over  the  sky.  Thus  the  space  occupied  by  these  13  sources  is  more 
likely  to  be  a  sphere  whose  diameter  is  smaller  than  the  equivalent. 
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width  of  the  galactic  plane  which  is  about  600  pc  (M.  Schmidt  1962b). 

The  volume  of  such  a  sphere  is  11  x  lO^  pc3.  if  this  contains  13 
sources,  then  the  number  density  is  one  radio  star  per  1.0  x  10?  pc3  or 
d  =  200  pc,  or  r  =  d/2  =  100  pc  as  the  average  distance  of  the  sun  from 
the  nearest  radio  star.  An  equally  rough  estimate  for  the  lower  limit 
to  the  distance  is  implied  by  the  lack  of  detectable  proper  motion,  less 
than  o'.'05/annum  for  30  48  and  30  286.  If  the  effect  of  standard  solar 
motion  only  is  considered,  then  the  minimum  distance  of  both  objects 
must  be  near  60  pc.  Any  additional  space  motion  of  the  stars  would 
probably  increase  this  minimum  distance  by  an  unknown  amount. 

Both  of  these  estimates  are  rough  and  should  only  be  taken  to  in¬ 
dicate  an  order  of  magnitude  for  the  distance.  In  the  calculations  of 
Section  IV  we  have  used  a  mean  distance  to  the  three  radio  stars  of 
100  pc. 


VII.  CONCLUSIONS 

The  three  radio  sources  identified  here  seem  to  be  the  first  mem¬ 
bers  of  an  entirely  new  class  of  objects.  The  physical  conditions 
giving  rise  to  the  radio  and  optical  radiation  are  almost  entirely  un¬ 
known.  The  parameters  computed  by  the  conventional  synchrotron  theory 
in  Section  IV  lead  to  strange  conditions  such  as  high  magnetic  fields 
and  short  replenishment  times  for  the  energy  stored  in  the  electrons. 

If  these  conditions  actually  apply,  we  have  "stellar"  objects  which  are 
capable  of  accelerating  high  energy  particles  in  short  times  ( T =  E/L  » 
years  to  a  few  102  years)  by  an  unknown  process.  If  there  are  as  many 
as  50  such  radio  stars  in  our  Galaxy,  these  objects  may  be  an  important 
source  for  cosmic  rays.  Work  now  in  progress  on  additional  identifica¬ 
tions  should  lead  to  a  significantly  larger  sample  for  study. 
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Figure  1.  The  two  color  diagram  for  various  types  of  objects. 

Line  C  is  for  unreddened  main  sequence  stars ,  line 
B  is  for  black  body  radiation,  and  line  A  is  for 
sources  whose  energy  distribution  is  of  the  form 
F( v)  pC  v"n.  Values  of  n  from  0.0  to  2.0  are  marked 
along  curve  C.  Colors  for  the  three  radio  stars  are 
plotted  aa  filled  circles.  Colors  for  a  few  old 
novae  are  shovn  as  crosses.  The  tics  on  curve  B  in¬ 
dicate  values  tabulated  in  Table  A3* 


Figure  2.  Finding  chart  for  3C  48  taken  from  a  10®  exposure 

vith  the  200-inch .  Local  photometric  standard  stars 
B  and  D  are  marked.  The  data  are  V  -  13.53,  B-V  - 
0.50,  U-B  -  0.00  for  star  Bj  V  -  14.54,  B-V  -  0.66, 
and  U-B  ■  0.05  for  star  D. 


Figure  4.  Finding  chart  for  3C  286 
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Figure  5 


The  radio  power  spectrua  of  3C  48  taken  from  Table  3- 
The  tvo  straight  lines  have  slopes  of  -O.67  and  -O.47. 
The  slope  of  -O.67  is  that  required  in  Figure  8  to  fit 
the  optical  data. 
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Figure  7.  The  radio  power  spectrum  for  3C  286  taken  fro*  Table  5 
The  line  has  a  slope  of  -0.28. 


Figure  8.  Comparison  of  the  radio  and  optical  power  levels  for 
3C  k8  with  the  theory  of  equation  (5).  The  optical 
data  are  from  Table  6. 
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Figure  9»  Detailed  plot  of  the  optical  data  for  3C  48  (Table  6) 
compared  vith  the  theoretical  curves  from  equation  (j) 
for  various  cut-off  frequencies.  The  critical  fre¬ 
quencies  at  the  high  energy  cut-off  are  given  on  the 
curves  in  units  of  lO^1*  c/o- 
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Figure  10.  Comparison  of  the  radio  and  optical  power  levels  for 
3C  196  with  the  theory  of  equation  (5).  The  optical 
data  are  from  Table  6. 


Figure  11.  Comparison  of  the  radio  and  optical  power  levels  for 
3C  286  with  the  theory  of  equation  (5).  The  optical 
data  are  from  Table  6. 


Table  1 


Photometric  Data  for  the  Three  Radio  Stare 


Object 

Date 

V 

B-V 

U-B 

Remarks 

3C  48 

Oct.  23/24,  I960 

16.06 

0.38 

-0.61 

Nov.  19/20,  1960 

16.02 

0.48 

-0.61 

60  inch 

Jan.  12/13,  1961 

16.11 

0.42 

-0.61 

Jan.  13/14,  1961 

16.13 

0.39 

-0.61 

Jan.  14/15,  1961 

16.02 

0.49 

-0.60 

Jan.  16A7,  1961 

16.13 

0.40 

-0.59 

Aug.  17/18,  1961 

16. 31' 

O.40 

-0.52 

Oct.  11/12,  1961 

16.333 

0.340 

-0.579 

Taken  15  min- 

Oot.  11/12,  1961 

16.289 

0.393 

-0.555 

uter  apart 

Dec.  4/5,  1961 

16.44 

0.35 

-0.57 

Obe.  by  Bean 

Dec.  5/6,  1961 

16.40 

.0.42 

-0.64 

Obs.  by  Bare 

17.79 

0.57 

-0.43 

3C  286 

June  2/3,  1962 

17.25 

0.26 

""  "1 

0.91 

Table  2 


Radio  and  Optical  Positions  for  Three  Radio  Stars 


Source 

1 

a 

Remarks 

3C  48 

01h34m50?3 

01  34  49.82 

0?6 

0.04 

Radio 

Optical 

3C  196 

08  09  59.6 

08  09  59.41 

0.4 

0.05 

♦48  22  06 
+48  22  07.9 

2.6 

0.5 

Radio 

Optical 

HE1 

13  28  50.8 

13  28  50.74 

+30  45  55 
♦30  45  59.5 

Iff! 

Optical 

Table  3 


e 

Radio  Tlux  Measurements  for  3C  43 


Observer 

v (Me/s) 

T(v)* 

Referenoes 

R3S  01.01 

81 

80 

1 

2C  133 

81 

102 

— 

2 

BSS  92 

100 

60 

20 

3 

3C 

159 

50 

22 

4 

HBH 

159 

60 

B 

5  . 

ERL 

178 

43 

6 

CIT(K) 

475 

33 

Wtm 

7 

CIT(K) 

710 

26.0 

2 

7 

Heesohen 

750 

21.6 

8 

CIT{MM) 

958 

20.4 

5 

9,10 

GIT(K) 

1420 

16.1 

2 

7 

Goldstein 

1423 

17 

5 

11 

cit(b) 

2841 

UBS 

5 

12 

Heesohen 

3000 

K! 1 

9 

8 

CIT(K) 

3200 

8.1 

_ 

11 

7 

*  The  flux  unit  is  10*26  W  m‘2(o/s)_1. 

The  referenoes  are  given  at  the  bottom  of  Table  5. 


Table  4 


Radio  flux  Measurements  for  3C  196 


Observer 

v  (Mo/s) 

F(v)* 

urn 

cs 

Remarks 

Whitfield 

38 

125 

56 

13 

Confused  reoord 

R3  08.01 

81 

100 

— 

1 

2C  724 

81 

124 

— 

2 

BSS  84 

100 

70 

20 

3 

30 

159 

66 

30 

4 

HBH 

159 

40 

50 

14 

ERL 

178 

66 

10 

6 

CIT(K) 

710 

25.2 

1.2 

7 

Heesohen 

750 

25.7 

10 

8 

CIT(MM) 

958 

19.4 

4 

9,10 

CIT(K) 

1420 

14.9 

3 

7 

Goldstein 

1423 

17 

8 

11 

CIT(B) 

28U 

7.5 

5 

12 

Heesohen 

3000 

10 

8 

CIT(K) 

3200 

WBa 

20 

7 

*The  flux  unit  is  lO”2^  W  *"2(o/s)"1. 

+  The  referenoes  are  given  at  the  bottom  of  Table  5. 


Table  5 


Radio  Flux  Measurements  for  3C  286 


Observer 

v(Me/s) 

F(v)* 

e7(* ) 

References 

2C  1120 

81 

50 

•• 

2 

3C 

159 

30 

23 

4 

ERL 

178 

22 

10 

6 

CIT(K) 

710 

19.1 

5 

7 

CIT(MM) 

958 

18.8 

7 

9,10 

CXT(K) 

1420 

16 

25 

7 

Goldstein 

1423 

17 

6 

11 

CIT(B) 

2841 

12 

6 

12 

*  The  flux  unit  is  10“^  W  m~Ho/a)~^- 
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Table  6 


Optical  Flux  for  the  Three  Radio  Stare 


Object 

v  xlO*1^ 
e/s 

F(vf 

€p($) 

Remarks 

30  !fi 

8.32 

.669 

1 

0  -  16.01 

6.77 

.958 

1 

B  -  16.60 

5.42 

1.29 

1 

V  -  16.20 

8.07 

.718 

1  " 

7.09 

.900 

1 

5.94 

1.19 

1 

Baum  Eight 

5.46 

1.35 

1 

Color  Data 

5.31 

1.46 

1 

* 

Reduced  to 

4.50 

1.94 

3 

V  -  16.20 

3.99 

1.78 

4 

3.53 

2.24 

12 

3.04 

2.26 

14  J 

30  196 

8.32 

.124 

2 

0  .  17.93 

6.77 

.195 

2 

B  -  18.36 

5.  a 

.302 

2 

V  -  17.79 

3C  286 

8.35 

U  -  16.60 

6.81 

B  -  17.51 

5.44 

V  -  17.25 

*  The  unite  for  F(  v  )  an  10"29  V  ■''2(o/e)"1. 
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) 


♦The  upper  cut-off  frequencies  used  are: 
3C  48  v2  =  7xl014  c/s 
3C  196  v2  =  4x  1014  c/s 
3C  286  v2  =10 x  1012  c/s 


Table  8 


Soaling  Factors  for  Other  Distances 


r(pcs) 

My-Myd) 

V(r) 

V(l) 

L(r) 

I»(l) 

ESI 

mm 

ET(r) 

%(l) 

T(r) 

r(i) 

1 

0 

1 

1 

i 

1 

n 

10 

-5 

103 

102 

0.52 

2.7jc102 

ISl. 

lit 

104 

0.27 

7.2xlQ4 

2*2 

103 

109 

106 

0.14 

2.0xl07 

20 

104 

-20 

1012 

108 

0.072 

5.2xl09 

52 

O 

fH 

-25 

1015 

0.037 

_ 

1.4X1012 

140 

Table  9 


Average  Values  of  Physical  Parameters  at  100  pcs 


V 

(cm3) 

H 

(gauss) 

1 

L 

(erg  sec"1) 

r 

(year) 

L 

(erg  sec-1) 

r 

(year) 

Halo 

1049 

5xl0‘3 

4X1043 

SxlO29 

3000 

6xl031 

| 

Core  S»l" 

2x10^ 

5x1  O'2 

7X1041 

9x1 029 

100 

5xl031 

Eg 

s-o!i 

2xl039 

4X10*1 

4X1040 

9X1029 

5 

SxlO31 

APPENDIX  A 


Theoretical  U-B,  B-V  Colors  for  Spectral 
Distributions  of  Arbitrary  Form 

In  many  applications  it  is  desirable  to  compute  the  colors  on  the 
U,  B,  V  system  of  radiant  sources  with  given  energy  distributions. 

This  can  be  done  once  the  transmission  functions  s(X)  of  the  U,  B,  V 
system  are  known.  Melbourne  showed  (i960)  that  the  S(X)  functions  tabu¬ 
lated  by  H.  L.  Johnson  (1955)  as  representing  the  system  will  not  pre¬ 
dict  exactly  the  observed  U-B,  B-V  colors  for  real  stars  of  known  energy 
distributions  unless  a  systematic  zero  point  correction  of  0?17  is  added 
to  the  computed  natural  u-b  magnitude  and  a  correction  of  -0.13  is  ap¬ 
plied  to  the  computed  natural  b-v  magnitude  before  applying  Johnson's 
empirically  determined  transformation  equation  (5)  of  Ap.  J.  117,  313, 
1953*  These  corrections  assume  that  there  is  no  color  equation  between 
the  theoretical  colors  based  on  S(X)  and  the  U-B  and  B-V  values  adopted 
for  standard  stars  in  the  sky.  In  this  appendix  we  confirm  Melbourne's 
zero  point  procedure  for  B-V  colors  by  showing  that  his  adopted  S(X) 
functions  predict  the  natural  colors  (b-v)0  outside  the  atmosphere  for 
seven  real  stars  such  that  the  regression  line  of  (b-v)0  on  the  known 
B-V  colors  has  a  slope  of  1.00.  However,  this  is  not  the  case  in  U-B 
where  the  adopted  S(X)  functions  predict  (u-b)D  colors  with  a  longer 
baseline  than  the  observed  U-B;  i.e.  a  color  equation  exists.  The 
equations  derived  in  this  appendix  can  be  used  to  convert  theoretical 
calculations  based  pn  the  adopted  S(X)  to  the  empirical  U,  B,  V  system. 
Predicted  U-B,  B-V  colors  for  any  arbitrary  flux  distribution  function 
F(X)  can  be  obtained  therefrom. 

Table  Al,  taken  from  Melbourne’s  thesis,  lists  the  transmission 
functions  for  the  theoretical  photometric  system  designated  by  u,  b,  v. 
As  anticipated  in  the  last  paragraph,  this  system,  aside  from  zero 
point  corrections,  is  close  to,  but  not  identical  with,  U,  B,  V.  The 
S(X)  functions  were  constructed  by  A.  D.  Code  from  Johnson's  tabulated 
values  for  zero  air  mass  by  applying  the  monochromatic  extinction  of 
the  atmosphere  for  one  and  two  air  masses  as  determined  from  spectrum 
scanner  observations  (Code,  as  tabulated  by  Melbourne  i960,  Table  2). 

We  have  tabulated  these  functions  because  they  are  needed  in  what  fol¬ 
lows  and  they  appear  nowhere  in  the  literature.  The  S(X)  for  one  air 
mass  was  tabulated  by  Arp  (1961)  in  his  discussion  of  U-B,  B-V  colors 
for  black  bodies,  but  S(X)0  and  S(X)g  are  needed  for  a  full  treatment 
of  the  problem. 

The  S(X)]_,  and  S(X)g  functions  of  Table  Al  were  used  to  operate  on 
the  relative  energy  distributions  F(X)  for  seven  real  stars.  The  F(X) 
data  were  obtained  from  Minnaert  (1953/  Table  l).  Code  (i960)  and 
Melbourne  (i960).  Code  and  Melbourne  tabulate  2.5  log  F(v)  where  F(v) 
is  the  energy  flux  per  unit  frequency  interval.  The  F(v)  values  were 
converted  to  relative  energy  flux  per  wavelength  interval  F(X)  by  multi¬ 
plying  F(v)  by-X-2.  The  seven  stars  chosen  were  either  relatively  free 
from  Fraunhofer  lines  or,  as  in  the  case  of  the  sun  (Minnaert  1953, 

Table  1),  the  continuum  level  was  already  adjusted  to  compensate  for  the 
effect  of  the  lines. 


Theoretical  colors  on  the  natural  system  of  S(A)  for  one  and  two 
air  masses  were  computed  from 


(u-b^  -  2.5  log 


L 

I 


o  Sb(l)j  F(l)  dl 
Su(X)i  F(X)  dX 


1  ■  1  or  2, 


(Al) 


with  a  similar  equation  for  (b-v)^  and  (b-v)2.  In  an  obvious  notation 
the  subscripts  refer  to  1  or  2  air  masses.  Results  for  the  seven  stars 
are  listed  in  Table  A2  where  the  observed  U-B,  B-V  colors  are  shown  in 
columns  2  and  3>  the  computed  natural  colors  (b-v)^,  (b-v)2  in  columns 
4  and  5,  (u-b)^,  (u-b)g  in  columns  6  and  7,  and  the  differences 
k12  s  (t>-v)2"(^*v)l  ana  ^56  *  (u-b )2- (u-b ) j.  in  columns  8  and  9-  Note 
that  k^2  a°d  keg  represent  the  color  extinction  between  sec  Z  »  2  and 
sec  Z  «  1  and  therefore  are  related  to  the  usual  extinction  coefficients 
defined  by 

k12  s  kl"k2  (*>-v)o  (A2) 


k5 6  s  k5-k6  •  (A3) 

As  was  expected,  the  extinction  coefficients  k12  and  *56  ar®  color  de¬ 
pendent  according  to  our  calculations  (l.e.  k2  and  kg  are  not  zero). 
From  Table  A2  it  is  easily  shown  that 


k12  =  0.102-0.024  (b-v)Q  (A4) 

which  is  in  excellent  agreement  with  observational  determinations  of 
extinction  on  good  nights  over  the  past  10  years  on  Mount  Wilson.  The 
extinction  coefficient  k^g  cannot  be  represented  by  such  a  simple  equa¬ 
tion  over  the  range  of  (u-b)0  covered  by  our  seven  calibrating  stars. 

The  value  of  kg  depends  on  (u-b)0  itself,  which  shows  that  equation  (A3) 
does  not  represent  the  broad  band  ultraviolet-blue  extinction.  The 
variation  of  kg  with  (u-b)0  has  been  found  observationally  (Johnson, 
unpublished  communication).  It  should  be  emphasized  that  in  setting  up 
the  U-B  system  in  1953>  Johnson  intentionally  put  k6  =  £  to  avoid  the  . 
great  complication  involved  in  the  reduction  procedure  of  k6  =  f[(u-b)gj. 
Therefore,  if  we  are  to  compare  theoretically  computed  u-b  indices  with 
U-B  values  assigned  to  real  stars,  we  must  ignore  the  color  variation 
of  k *6  even  though  k6  +  0  in  nature.  The  adopted  U-B  values  for  stars 
in  the  sky  therefore  contain  slight  inconsistencies  with  a  purely  theo¬ 
retical  system.  Consequently,  in  the  following  we  shall  relate  the 
computed  (u-b)i  values  directly  to  U-B  rather  than  use  the  theoretical 
keg  value  from  Table  A2  to  obtain  (u-b)0  from  (u-b )i  in  the  usual  way 
[i.e.  (u-b)0  =  (u-b)j^  -  k  =  6j  and  then  to  find  (u-b)0  =  f(U-B).  We  tried 
this  latter  procedure,  which  is  theoretically  correct  because  it  takes 
the  non-zero  kg  term  into  account,  but  the  correlation  of  (u-b)Q  and  U-B 
was  non-linear,  deviating  by  ±  0?05  from  straight  line.  The  foregoing 
explanation  is  believed  to  be  correct  (i.e.  small,  residual  inconsis¬ 
tencies  in  the -actual  U-B  values),  because  the  non-linearity  was  elimi¬ 
nated  in  a  plot  of  (u-b)i  vs  U-B.  Conversations  with  Harold  Johnson 
after  these  calculations  were  completed  confirm  that  this  is  undoubtedly 
the  case. 


With  the  foregoing  precepts  and  with  the  data  of  Table  A2,  the 
following  extinction  and  color  equations  can  be  derived.  From  equation 
(A4)  and  from  the  definition  of  extinction  coefficients,  which  ie 

(b-v)0  =  (b-v)8ec  z  -  (b-v)0J  sec  Z,  (A5) 


it  follows  that 


(b-v) 0  =  1.024  (b-v)!  -  0.104  .  (A 6) 

• 

With  (b-v)0  now  known  by  using  equation  (A6)  and  the  (b-v)^  values 
of  Table  A2,  we  correlate  (b-v)0  with  the  observed  B-V  for  the  seven 
calibrating  stars.  The  result  is 

B-V  =  1.00  (b-v)0  +  0.91  (A7) 

or 

B-V  =  1.024(b-v)1  +  0.8l  (a8) 

which  shows  that  there  is  no  color  equation  between  Table  A1  and  the 
B-V  system  if  the  extinction  is  treated  in  this  way. 

The  U-B  situation  is  different.  No  extinction  values  were  applied 
to  (u-b)i  to  obtain  (u-b)0  for  the  reasons  previously  discussed.  Direct 
comparison  of  (u-b)]^  and  U-B  from  Table  A2  gives  the  well  determined  re¬ 
lation 

U-B  =  0.921  (u-b)!  -  1.308  (A9) 

which  is  the  color  equation  between  the  observed  U-B  system  and  the 
theoretical  (u-b)^  system. 

An  objection  might  be  raised  to  our  procedure  in  that,  in  actual 
fact,  we  may  not  know  the  true  F(A)  for  the  seven  calibrating  stars  of 
Table  A2  due  to  the  effects  of  line  blanketing,  with  the  result  that 
our  predicted  (u-b)^  colors  may  be  in  error.  While  this  objection  may 
be  true  for cr Boo,  and  to  a  lesser  extent  for  a  Lyrae  with  its  strong 
hydrogen  lines,  it  is  not  true  for  the  other  stars.  10  Lac  and  ^UMa 
are  so  hot  as  to  have  few  lineBj  HD  19^5  and  HD  140283  are  cool  sub¬ 
dwarfs  with  extremely  weak  Fraunhofer  lines  and  the  blanketing  effect 
is  almost  nil,  while  the  effect  in  the  sun  has  been  corrected  for  by 
Minnaert  before  giving  his  F(X)  of  his  Table  1  (1953).  Therefore  we 
interpret  the  slope  coefficient  in  equation  (A9)  to  differ  significant¬ 
ly  from  1.00  and  suggest  that  a  color  equation  exists.  We  are  not  com¬ 
pletely  confident  of  the  result  because  of  the  blanketing  problem,  but 
accept  it  for  the  present. 

Equations  (A8)  and  (A9)  solve  the  problem.  To  compute  U-B  and  B-V 
for  any  arbitrary  energy  distribution  F(l),  proceed  as  follows. 

(l)  Use  the  S(l)i  functions  of  Table  A1  and  use  equation  (Al)  to  com¬ 
pute  natural  colors  (u-b)^  and  (b-v)^.  (2)  Use  equations  (A8)  and  (A9) 

to  find  U-3  and  B-V.  A  check  on  the  internal  consistency  of  this  pro¬ 
cedure  is  made  by  applying  the  equations  to  the  seven  calibrating  stars. 
Comparison  of  the  computed  colors  in  columns  (10)  and  (ll)  of  Table  A2 
with  the  observed  colors  of  columns  (2)  and  (3)  shows  good  agreement. 


The  above  procedure  differs  from  that  used  by  Arp  (1961)  in  that 
we  have  used  a  color  equation  (A9)  to  the  U-B  system  whereas  Arp  had 
none.  Our  zero  points  also  differ  slightly.  Table  A3  shows  the  dif¬ 
ference  between  Arp’s  calculations  of  the  black  body  line  and  our  own. 
The  difference  is  very  small. 

In  the  theory  of  synchrotron  radiation,  energy  distributions  of  the 
form  F(v)e<v-n  (flux  per  unit  frequency  interval)  are  found.  Calcula¬ 
tion  of  U-B  and  B-V  values  have  been  made  for  these  distributions  for 
various  values  of  n  ranging  from  0.0  to  2.0  in  steps  of  0.2.  The 
sources  were  assumed  to  radiate  as  FfA)-*!11"2  (per  unit  wave  length 
interval).  The  results  are  given  In  Table  A  4  and  have  been  discussed 
in  the  body  of  the  paper. 


TABLE  A1 


Adopted  Transmission  Functions  After  Melbourne  and  Code 


■Hi 

S(A)0 

S(A)i, 

One  Air  Mass 

S(X)2,  Two  Air  Masses 

bo 

Vo 

Ui 

bi 

Vi 

ua 

ha 

wtm 

30 

0.025 

9 

31 

0.250 

0.060 

32 

0.680 

IJBIB  r  »■ 

33 

1. 137 

Bw  r 

34 

1.650 

Biw.  r 

35 

0.000 

2.006 

0.000 

1.000 

36 

0.015 

2.250 

0.006 

1.250 

0.000 

37 

0. 100 

2.337 

0.080 

0.  040 

38 

0.500 

1.925 

0.337 

0.  250 

39 

.  1.800 

0.  650 

1.425 

0.600 

0.870 

40 

3.620 

0.  197 

2.253 

0.140 

1.745 

41 

3.910 

0.070 

2.  806 

0.030 

2.025 

42 

4.  000 

0.000 

2.950 

0.000 

2.200 

43 

3.980 

3.  000 

2.290 

44 

3.780 

2.937 

2.290 

45 

3.500 

2.780 

2.200 

46 

3. 150 

2.520 

2.050 

47 

2.700 

0.000 

2.230 

1.840 

0.000 

48 

2.320 

0.020 

1.881 

0.020 

1.580 

0.020 

49 

1.890 

0.280 

1.550 

0. 175 

1.325 

0.  175 

50 

1.530 

1.180 

1.275 

0.900 

1.079 

0.830 

51 

1.  140 

2.170 

0.975 

1.880 

0.840 

1.570 

52 

0.750 

2.970 

0.695 

2.512 

0.580 

2.150 

53 

0.500 

3.300 

0.  430 

2.850 

0.360 

2.420 

54 

0.250 

3.250 

0.210 

2.820 

0. 180 

2.  425 

55 

3.000 

0.055 

2.625 

0.050 

2.270 

56 

2.700 

0.000 

2.370 

0.000 

2.035 

57 

2.320 

2.050 

1.780 

58 

2.000 

1.720 

1.530 

59 

1.670 

1.413 

1.255 

60 

1.280 

• 

1.068 

0.985 

61 

0.920 

0.795 

0.745 

62 

0.650 

0.567 

0.505 

63 

0.430 

0.387 

K  IHil 

64 

0.260 

0.250 

■i  wyXwr. 

65 

0. 175 

0. 160 

■i'  itlm 

66 

0.125 

0.110 

0. 100 

67 

0.100 

0.081 

0.080 

68 

0.070 

0.061 

0.055 

69 

0.050 

0.045 

0.030 

70 

0.030 

0.028 

0.025 

71 

0.020 

0.017 

0.010 

72 

0.010 

0.007 

0.000 

Tabulated  are  the  transmissions  of  two  aluminum  reflections,  filters,  1P2T 
photomultiplier,  and  atmosphere  for  zero,  one,  and  two  air  masses. 


Table  A2 


Computed  and  Observed  Colore  for  the  Seven  Calibrating  Stare 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

Star 

(®-V)obe 

<D-«obe 

(b-v^ 

ma 

(«). 

10  Lao 

-0.20 

-1.04 

-0.994 

pya 

-1.05 

77  0  Ma 

-0.21 

-0,68 

-0.984 

-0.857 

0.683 

1.121 

0.127 

Mr*?! 

-0.20 

-0.68 

HD  19445 

O.46 

-0.24 

-0.338 

-0.223 

1.182 

1.596 

0.115 

0.414 

O.46 

-0.22 

HD  140283 

0.48 

-0.20 

-0.288 

-0.174 

1.161 

1.567 

0.114 

0.406 

0.51 

-0.24 

a  Boo 

0.37 

-0.09 

-0.436 

-0.320 

1.332 

1.725 

0.116 

0.393 

0.36 

-0.08 

a  Lyr 

0.00 

-0.01 

-0.839 

-0.714 

1.410 

1.762 

0.125 

0.352 

-0.05 

-0.01 

Sun 

0.62 

0.10 

-0.196 

-0.090 

1.512 

1.933 

.0.106 

0.421 

0.61 

0.08 

Table  A3 


Black  Body  Line  in  the  U-B,  B-V  Diagram 


■ 

Present 

Arp 

1  1 

tJ-B 

B-V 

■rejf 

msm 

CO 

-1.28 

-0.44 

-1.33 

-0.46 

25,000 

-1.13 

-0.21 

-1.17 

-0.23 

20,000 

-1.06 

-0.15 

-1.09 

-0.17 

12,000 

-0.83 

+0.05 

-0.84 

+0.04 

8,000 

-0.53 

+0.35 

-0.52 

+0.34 

6,000 

-0.26 

+0.63 

-0.22 

+0.61 

5,000 

-0.10 

+0.79 

-0.05 

+0.78 

4,000 

+0.36 

♦1.13 

♦0.37 

+1.12 

3,300 

♦I.44 

+0.83 

+1.44 

3,000 

+1.67 

+1.14 

_ 

+1.66 

APPENDIX  B 


Conversion  of  UBV  Optical  Data  to  Absolute  Flux  Units 

The  most  recent  calibration  of  the  energy  flux  in  absolute  units 
arriving  from  a  star  of  apparent  visual  magnitude  V  «*  0  is  due  to 
Willgtrop  (i960).  We  adopt  Willstrop's  Table  1  for  the  energy  flux 
per  A  at  A5390  for  stars  of  V  =  0.00,  and  take  the  known  relative 
#energy  distribution  curves  for  stars  of  different  B-V  values  to  compute 
the  absolute  flux  calibrations  at  the  effective  wavelengths  of  the  U, 

B,  and  V  points  for  stars  of  U  =  B  =  V  =  0. 

a )  Effective  Wavelengths  for  the  UBV  Filters 


The  effective  wavelength  of  any  filter  system  depends  upon  the 
energy  distribution  F(A)  of  the  source  which  is  measured  with  the 
system.  We  define  the  effective  wavelength  in  the  usual  way  by 


A  = 


r 


AS (A)  F(A)  dA 
F(A)  dX 


r  S ( A) 


(Bl) 


and  compute  A  for  sources  characterized  by  F(A)ocAn-®. 


The  effective  wavelengths  calculated  from  equation  (Bl)  using  the 
S ( A )  function  from  Table  A1  are  for  the  u,  b,  v  natural  'system,  which 
1  s  related  to  the  UBV  system  via  equations  (a8 ) ,  and  (A9)  of 

Appendix  A.  From  these  equations  and  from  Au,  A^,,  Av  for  one  air  mass, 
we  can  obtain  Ay,  XB,  Xy  for  the  actual  UBV  system  as  follows. 


The  relation  between  two  color  systems,  characterized  by  effective 
wavelengths  Ax,  Ag,  and  A3,  X4,  when  the  measurements  are  made  on  black 
bodies  is 


Ccr)v, 


1 _ L 

Ai  A*. 

X7T 

As  A* 


(B2) 


For  sources  ve  are  considering  where  F(A)ocxn-2 
will  be 


> 


the 


correct  relation 


(B3) 


but  for  color  systems  where  the  effective  wavelengths  are  very  close  to 
one  another,  equation  (B2)  and  ( B3 )  give  nearly  identical  answers.  We 
use  equation  (B2)  in  the  following. 

The  effective  wavelength  of  V  wais  obtained  by  assuming  Ay  a  Xv  for 
SV(A)0  in  equation  (Bl).  This  gives  Ay  listed  in  column  6  of  Table  Bl 
for  various  n  values.  The  Xg  follows  from  equation  (B2)  and  equation 


(A8)  as 


where  X^  is  the  effective  wavelength  of  the  natural  b  system  for  1  air 
mass  and  likewise  for  v  in  Xvj_.  A  cheek  on  the  above  procedure  is 
available  from  equation  (A7)  which  requires  that  Xg  =  XgQ.  The  check 
was  made  and  was  found  to  be  excellent.  Column  5  of  Table  B1  gives  Xg. 


For  the  ultraviolet,  equation  (A9)  requires 


_1_ 

Xy 


(B5) 


where  all  quantities  on  the  right  side  are  now  known  or  were  computed 
from  equation  (Bl).  Column  4  of  Table  B1  gives  Xy  for  the  various  n 
values. 


b )  Calibration  of  Magnitudes  in  Absolute  Flux  Units 


Willstrop,  in  his  Table  1  (i960),  gives  absolute  monochromatic 
fluxes  at  X  «  5390  S,  a  wavelength  which  differs  from  Xy.  We  need  the 
ratio  F(Xv)/f(5390)  to  correct  Willstrop’s  values  to  monochromatic 
fluxes  at  Xy.  For  the  sources  we  are  considering,  where  F(X)  =  CXn_  , 
it  follows  that 


log 


F(  Xy  ) 
F(5390) 


(n-2) 


(B6) 


From  equation  (b6),  from  Willstrop' s  Table  1.  and  from  the  conver¬ 
sion  of  flux  units  in  ergs  cm-2  sec'^  A*1  to  W  m'2(c/s)_1  by  multiplying 
the  former  by  105  X2/c  to  get  the  latter,  we  obtain  the  flux  at  Xy  for  a 
source  of  visual  magnitude  V  =  0  tabulated  in  column  9  of  Table  Bl, 
which  i 8  one  of  the  three  answers  we  require.  The  flux  at  the  effective 
wavelengths  of  the  B  and  U  bands  are  now  found  as  follows.  For  a  given 
index  n  (for  example  n  =>  2),  the  flux  in  the  source  at  the  wavelengths 
Xy  and  XB  are  in  the  ratio  (vy/^’2  =1.5596  (for  n  =  2)  per  unit  fre¬ 
quency  interval.  Thus,  if  the  visual  magnitude  of  the  source  were  0.00, 
the  absolute  flux  density_in  W  m*2  (c/s)"l  at  Xg  =  4440  A  would  be 
1.5596  times  the  flux  at  Xy  which  we  have  previously  calculated  and 
listed  in  Table  Bl.  But  Table  A4  or  Bl  show  that  B-V  =  +O.58  for  n  =  2, 
or  B  =  +O.58  by  the  requirement  that  V  ■  0.00.  Hence,  the  flux  at 
XB  s  4440  for  a  source  of  index  n  »  2  corresponding  to  a  blue  magnitude 
0.00  will  be  1.84  x  10‘23  W  m-2  (c/s)  .  Similar  calculations  were  made 

for  all  cases  and  are  tabulated  in  Table  Bl  and  these  numbers  solve  the 
problem. 


Table  B1 


Adopted  Effective  Wavelengths  and  Conversions  to  Absolute 
Flux  Units  for  the  U,  B,  V  System  for  Souroes  which 
Radiate  as  ?{V  )  »  CV"n 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

n 

B-V 

\ 

El 

FU)0 

0.0 

-0.91 

0.10 

3585 

4383 

5502 

1.79 

4.13 

3.81 

0.2 

-0.86 

0.15 

3586 

4387 

5505 

1.80 

4.15 

3.82 

0.6 

-0.77 

0.25 

3593 

4408 

5515 

1.81 

4.18 

3.84 

1.0 

-0.68 

0.34 

3598 

4412 

5525 

1.82 

4.21 

3.87 

1.4 

-0.60 

O.44 

3604 

4431 

5534 

1.83 

4.26 

3.90 

2.0 

-O.46 

0.58 

3610 

4440 

5545 

1.84 

4.31 

3.94 

U-B  and  B-V  taken  from  Table. A4  of  Appendix  A. 

The  unit  on  the  flux  is  10**’  W  m"*(c/s)-1  for  a  star  of  U»^V»0 
apparent  magnitude. 
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ADDENPUM  TO  THIS  PAPER 
VIII.  3C  48  AS  A  GALAXY 


After  this  paper  had  been  submitted,  new  evidence  strongly  sug¬ 
gests  that  3C  273 >  another  similar  object,  has  an  appreciable  redshlft 
(2  =  O.I58,  Schmidt  1963) .  The  spectrum  of  3C  48  can  also  be  under¬ 
stood  if  it  also  has  a  redshlft  (Greenstein  and  Matthews  1963).  This 
section  has  been  included  to  discuss  the  observational  data  for  3C  48 
under  the  assumption  that  it  is  a  galaxy. 

A.  The  redshlft  of  1  »  O.3675  for  3C  48  (Greenstein  and  Matthews  19&3) 
can  be  interpreted  as  a  distance  of  1100  Mpc  using  a  Hubble  constant  of 
H  ■  100  km/sec  (Mpc)"^.  Using  this  distance  to  calculate  the  absolute 
magnitude  of  the  optical  object  gives  the  value  in  Table  9  which  has 
been  corrected  for  the  effects  of  the  redshlft*  but  not  for  interstel¬ 
lar  absorptions,  except  for  Cygnus  A.  The  absolute  magnitudes  for 
3C  48  and  3C  273  are  compared  to  other  more  normal  radio  galaxies  3C  295 
and  Cygnus  A  and  to  the  brightest  field  and  cluster  galaxies  known 
(Humason  et  al  1956).  3C  48  is  about  2  magnitudes  brighter  optically 
than  the  brightest  field  and  cluster  galaxies,  and  2.6  magnitudes 
brighter  than  the  average  of  3C  295  and  Cygnus  A.  Table  9  shows  that 
3C  273  is  even  brighter. 

The  observed  radio  and  optical  size  of  3C  48  of  4.  1  second  of  arc 
means  that  the  radio  and  optical  flux  is  coming  from  a  diameter  of 
45500  pc.  In  the  case  of  3C  273  (Hazard  et  al  1963;  Schmidt  1963), 
the  diameter  of  the  core  of  component  B  is  4  1000  pc.  The  faint  nebu¬ 
losity  visible  around  3C  48  is  66  kpc  long  if  at  1100  Mpc.  The  nebu¬ 
losity  is  probably  not  a  faint  underlying  galaxy  for  several  reasons. 

(l)  It  1 s  much  larger  than  one  would  expect  for  a  galaxy  at  that  dis¬ 
tance.  (2)  The  appearance  of  the  nebulosity  does  not  resemble  a 
galaxy  in  outline  or  in  the  distribution  of  brightness.  (3)  The  bright 
object  having  a  stellar  appearance  on  direct  plates  is  not  centrally 


*The  bolometric  magnitude  must  be  corrected  for  the  K  effect  due  to  the 
shift  of  the  spectrum  under  the  measuring  bands.  This  is  composed  of 
two  terms,  one  due  to  the  change  of  band  width  of  the  filters  relative 
to  the  spectrum  being  measured,  and  the  other  due  to  the  wavelength  de¬ 
pendent  effect  of  reducing  hertrochromatie  magnitudes  to  bolometric 
(see  eq.  B7  of  Humason  et  al  1956).  For  an  effective  index  of  n  -  2, 
the  second  term  is  zero.  The  first  term  is  always  2.5  log  (l+z). 

Hence  since  n  a  1.5  for  3C  48,  we  neglect  the  OV34  for  3C  48.  For 
3C  295  and  Cygnus  A,  the  K  correction  of  Humason  et  al  Table  Bill  was 
applied.  These  corrections  are  only  approximate  for  two  reasons: 

1)  only  the  first  term  (cz/h)  in  the  distance  equation  (Sandage  1961, 
eq.  50)  has  been  used;  and  (2)  the  selective  K  term  is  not  known  exactly. 
However,  they  suffice  for  our  present  purpose.  Except  for  Cygnus  A,  no 
correction  for  absorption  in  our  ovn  galaxy  has  been  applied. 
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located  (see  Section  II).  It  is  16  kpe  from  the  center  of  the  visible 
nebulosity.  (4)  The  integrated  apparent  magnitude  of  the  nebulosity 
omitting  the  "star"  is  mv  G4  16.7.  Its  absolute  magnitude,  given  in 
Table  9,  is  brighter  than  the  brightest  of  known  cluster  galaxies,  and 
1  magnitude  brighter  than  the  average  of  3C  295  and  Cygnus  A  which  are 
normal  radio  galaxies.  In  the  case  of  more  normal  and  radio  galaxies, 
the  nucleus  of  the  galaxy  is  usually  the  brightest  part  and  contributes 
significantly  to  the  absolute  magnitude.  The  integrated  absolute  mag¬ 
nitude  for  the  nebulosity  around  3C  48  does  not  include  the  effect  of 
any  such  nucleus. 

If  the  nebulosity  was  produced  by  some  phenomenon  (say  an  explo¬ 
sion)  at  some  time  in  the  past,  then  from  its  visible  extent  we  can 
infer  that  the  event  occurred  ^  1.8  x  105  years  ago.  The  minimum  cor¬ 
responds  to  assuming  that  the  expansion  has  occurred  at  the  velocity  of 
light  and  there  are  no  projection  effects. 

B.  Physical  Parameters.  The  physical  parameters  for  3C  48  at  a  dis¬ 

tance  of  1.1  x  109  pc  have  been  computed  in  exactly  the  same  manner  as 
in  Section  IV.  The  results  are  given  in  Table  10.  Note  that  the  value 
of  the  magnetic  field  for  minimum  total  energy  is  near  10-3  gauss,  much 
lower  than  suggested  in  Section  IV,  and  not  far  from  other  radio 
galaxies  (Maltby  et  al  1963).  The  total  energies  required  are  also 
similar  to,  but  lower  than,  some  of  the  other  radio  galaxies.  However, 
note  that  the  time  scales  are  all  very  short  compared  to  the  1.8  x  10? 
years  suggested  above.  We  can  resolve  this  difficulty  in  one  of  two 
ways.  (l)  The  nuclear  region  of  the  galaxy  has  been  continuously  pro¬ 
ducing  high  energy  particles  throughout  its  life  time.  The  total  energy 
required  since  birth  is  then  5  x  1059  erg  (v2  =  1010)  or  3  x  IO61  erg 
(v2  =  5  x  10*M .  (2)  The  magnetic  field  is  lower  by  at  least  the 

amount  required  to  make  the  lifetime  equal  to  1.8  x  105  years.  The 
physical  parameters  required,  under  this  latter  hypothesis,  are  given 

in  the  last  row  of  Table  10.  There  is  no  dependence  on  the  size  of  the 
object.  The  total  energy  required  in  both  these  cases  is  becoming  un¬ 
comfortably  large  and  it  will  be  even  larger  if  H  is  made  any  smaller. 

• 

C.  Optical  Fluctuations.  The  observation  of  the  optical  flux  varia¬ 
tions  (Section  III)  becomes  more  important  if  these  objects  are  galaxies, 
since  we  normally  do  not  thick  of  the  light  from  galaxies  as  varying.  A 
time  scale  of  6  months  for  the  variations  implies  a  size  of  the  emitting 
region  of  0.15  pc  on  the  basis  of  light  travel  time  across  the  object. 
If  we  are  seeing  statistical  variations  of  a  number  of  similar  events, 
about  10  such  events  would  be  happening  simultaneously  and  the  emission 
region  would  be  at  least  10  times  larger.  The  variations  could  still  be 
due  to  decay  of  optical  synchrotron  electrons  in  a  magnetic  field  of  ~  1 
gauss  (see  Section  III).  This  is  not  inconsistent  when  we  note  that 
the  total  energy  depends  very  little  on  which  hypothesis  we  use  in  the 
previous  subsection.  Thus  without  changing  the  total  energy  require¬ 
ments  very  much,  we  can  choose  the  magnetic  field  and  adjust  the  emit¬ 
ting  volume  to  fit.  The  discussion  of  the  lack  of  radio  fluctuations 
would  be  the  same  as  in  Section  III. 


Greenatein , 
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TABLE  9 


Absolut*  Magnitudes  Based  on  E  *  100  kny/seo  10'  pae 


Object 

OS 

(lao/aeo) 

”V 

< 

3C  48 

110,000 

16.2 

-24.3 

Neb  3C  48 

110,000 

18.7 

-22.8 

30  273 

47,600 

12.6 

-25.6 

Cyg  A 

16,830 

14.1 

-22.1 

3C  295 

138,000 

20.9 

-21.3 

Field  and 
Cluster1' 

— 

— 

-22.3 

- 

^  On  the  system  of  the  first  brightest. 

^Not  inoluding  K  correction.  For  Cyg  A  includes 
1?S6  absorption  and  0®76  incompleteness. 

*  Inoludea  K  correction. 


